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We study the breakup of HJ exposed to super-intense, femtosecond laser pulses with frequen¬ 
cies greater than that corresponding to the ionization potential. By solving the time-dependent 
Schrodinger equation in an extensive field parameter range, it is revealed that highly nonresonant 
dissociation channels can dominate over ionization. By considering field-dressed Born-Oppenheimer 
potential energy curves in the reference frame following a free electron in the field, we propose a 
simple physical model that characterizes this dissociation mechanism. The model is used to predict 
control of vibrational excitation, magnitude of the dissociation yields, and nuclear kinetic energy 
release spectra. Finally, the joint energy spectrum for the ionization process illustrates the energy 
sharing between the electron and the nuclei and the correlation between ionization and dissociation 
processes. 


PACS numbers: 33.80.Rv, 33.80.Gj, 82.50.Kx 

Molecular breakup processes induced by strong 
light-matter interactions are of fundamental inter¬ 
est. Molecular-specific ionization phenomena include 
charge-resonance-enhanced ionization [1], subcycle mul¬ 
tiple ionization bursts [2], electron-nuclei energy-sharing 
in above-threshold dissociative ionization [3], above¬ 
threshold Coulomb explosion [4], and interplay between 
multiphoton and tunneling ionization [5]. Dissocia¬ 
tion phenomena include above-threshold dissociation [6], 
bond-softening [7], bond-hardening [8, 9], and rescatter¬ 
ing induced dissociation [10]. All these processes were 
discovered in the low-frequency regime where absorption 
of several photons is needed to reach a breakup channel. 
Depending on the laser parameters, either dissociation or 
ionization dominates [11-13]. 

With advancements in light-source technology, 
extreme-ultraviolet (XUV) laser pulses of femto- and 
subfemtosecond duration are now produced from high- 
order harmonic [14-17] or free-electron lasers [18, 19]. 
New focusing techniques [20-22] led to XUV femtosec¬ 
ond (fs) pulses with peak intensities / > 4 x 10^^ W/cm^ 
[23]. Intense XUV pulses were, e.g., applied on rare 
gases to study sequential versus non-sequential multiple 
ionization [24, 25], and creation of charge states up to 
21 in Xe was observed [21]. For molecules, experiments 
on HeH^ [26] and N 2 [27] provided benchmark data 
for theory. The high-frequency regime is defined in 
this work as the regime where one-photon ionization 
is allowed. Since the photon resonance is much closer 
to the threshold for ionization than for dissociation, 
ionization is at first glance expected to dominate. 

Here we characterize some hitherto unobserved molec¬ 
ular breakup phenomena in the regime of super-intense, 
high-frequency, fs pulses, supplementing the phenomena 
known from the low-frequency regime, and adding in¬ 
sight to the general field of strong laser-matter interac¬ 
tion. The characteristics we Hnd include (a) even with 


full inclusion of nuclear dynamics, stabilization against 
ionization occurs, i.e., the ionization yield does not nec¬ 
essarily increase with intensity [28, 29]; (b) a mechanism 
by which dissociation, in contrast to the expectation from 
energy considerations, completely dominates over ioniza¬ 
tion; (c) control over the vibrational distribution, disso¬ 
ciation yield, and nuclear kinetic energy release (NKER) 
spectra by the parameters of the laser pulse; (d) insight 
into the energy sharing between electronic and nuclear 
degrees of freedom, as displayed by the joint energy spec¬ 
trum (JES). 

We consider HJ, a molecule that was used to discover 
many low-frequency processes that were later observed in 
more complex molecules (see, e.g., [3, 30]). The possible 
extension of our findings to other molecules is discussed 
in the end. We consider a co-linear model that includes 
the dimension that is aligned with the linearly polar¬ 
ized pulse. This model allows us to capture the essential 
physics, and to perform calculations over an extensive 
laser parameter regime. The time-dependent Schrodinger 
equation (TDSE) reads (atomic units are used through¬ 
out), 

i?,t) = [Te + T^+ VeNix,R)+ Vi{t)]'ii{x,R,t) 

( 1 ) 

with X the electronic coordinate measured 
with respect to the center of mass of the 
nuclei, R the internuclear distance, Tg = 
-(1/2^)979x2, Tn = -(l/mp)92/9i?2, Vg^ix,R) = 
-l/7(x - 77/2)2 -H aiR)-l/y/{x + R/2)^ + aiR) + l/R, 
and Vi{t) = —i(3A{t)d/dx, with rup the proton mass, 
/i = 2m-pj(flm-p -I- 1), /? = (™p + l)/''7rp, and a[R) 
the softening parameter producing the exact IsUg 
Born-Oppenheimer (BO) curve [3]. The initial ground 
state has the energy Eg = —0.5973, the dissociation 
limit Ed = —0.5, the equilibrium internuclear distance 
i?o = 2.064, and the ionization potential Ip = 1.1 at 
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FIG. 1. (color online). Continuum probabilities for inter¬ 
acting with intense laser pulses containing 100 cycles, corre¬ 
sponding to r = 2.1 — 15.8 fs. (a) and (b) DI and dissociation 
probabilities for lj = 2.278, (c) and (d) DI and dissociation 
probabilities for different and Fq. The dashed line in (c) 
traces the position of the largest DI rate, calculated from 
HFFT (see text). The dashed, straight line in (d) corresponds 
to oo = F^juP' — 2.41. The inset in (d) shows the nuclear dy¬ 
namics for pulse parameters tracing this line (see Fig. 2 and 
the accompanying discussion). 


i?o- We use vector potentials A{t) = (Fq/uj) g{t) cosojt, 
with Fq related to the peak intensity by / = Fq. 
The envelope is g{t) = exp (—41n(2)t^/T^), with r 
the FWHM, and the number of cycles defined by 
T = dTrVin 2Nc/ uj. All pulses used satisfy the non- 
relativistic criteria <C 1 [31, 32] and the dipole 

condition Up/2ujc <C 1 [33], with Up = Fq/Alu'^ the 
ponderomotive potential. 

We solve (1) numerically [34] and obtain yields for the 
dissociation and dissociative ionization (DI) channels of 
(see Suppl. Material [35]). Fig. 1(a) shows the DI 
probability Pion of resulting from laser pulses with 
T = 11.1 fs and oj = 2.278. Pion increases with Fq for 
lower amplitudes until it reaches a maximum at Fq ~ 6.2, 
whereafter it decreases. The latter behavior indicates 
stabilization with respect to DI [28]. Figure 1(b) presents 
the dissociation probability Pdis- For the lower values of 
Fq, we observe no dissociation as expected from pertu- 
bation theory. At Fq ~ 9, dissociation sets in, increasing 
with Fq, until it “saturates” at Fq ~ 15. 

To obtain a complete picture. Figs. 1(c) and 1(d) 
present Pjon and Pdis, as functions of Fq and For fixed 
UJ the suppression of DI for large Fq is evident [Fig. 1(c)]. 
Lobes corresponding to large Pjon are seen emanating 
from the origin. By using a version of high-frequency Flo- 
quet theory (HFFT) [36] generalized to include nuclear 
motion [35], we find that the largest lobe in Fig. 1(c) is 
along the largest peak in the DI rates. The dashed line 
traces the largest HFFT rate, in reasonable agreement 
with the TDSE result. The oscillatory behavior in Pion 
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FIG. 2. (color online). Schematic of the nuclear dynamics in 
the KH frame for different oo and r. For (a)-(c), ao = 3.25, 
and for (d)-(f), oq = 1.87. In each panel, from bottom to top, 
the lowest field-dressed BO curve is shown for different times 
corresponding to ao(t) — 0, 1.87, (and 3.25). The lowest three 
field-dressed vibrational levels are indicated. The horizontal 
lines indicate Ed- The arrows sketch the pathways of the 
vibrational WP during the turn-on and -off of the pulse. The 
final position of the WP after the pulse is shown in dark blue. 

was discussed for reduced-dimensionality atomic systems 
[28, 37, 38], and was attributed to the two-center char¬ 
acter of dressed the atomic wave function. In the case of 
Pdis in Fig. 1(d), we observe appreciable dissociation for 
pulses with ao = Pq/w^ > 2.41 indicated by the dashed 
line. In this regime, Pion + Pdis = 1, implying unity prob¬ 
ability for breakup. The remnants of the lobes from DI 
are present in Fig. 1(d), an indication of the different time 
scales for DI and dissociation processes: DI occurs during 
the pulse, and what is not ionized, dissociates after the 
pulse. The physics of the dashed line will be explained 
below. 

In order to elucidate the origin of the onset of dissoci¬ 
ation in Fig. 1, it is instructive to transform Eq. (1) into 
the Kramers-Henneberger (KH) frame [39, 40], where the 
complete laser-molecule interaction is contained in the 
modified electron-nuclei interaction 14 n(;c -I- (3a(t),R), 
with the quiver motion a{t) = J* A{t')dt'. In the case 
of a monochromatic laser field, the Floquet ansatz for 
the wave packet (WP) and the Fourier expansion of 
KN(a; -b f3a{t),R) results in a coupled set of equations. 
For large uj, effectively only the zeroth order Fourier com¬ 
ponent, VQ{x,R,aQ) = (w/27r) VeN(x + /3a(t), R)dt, 
remains, resulting in the structure equation [28, 35, 40] 

[ffe(x, P; ao) + Tn] u(x, R; ao) = W (ao)u(a:, R; ao), 

( 2 ) 

with F[e{x, R; oq) = Pg + Vb(a;, P;ao) the field-dressed 
electronic Hamiltonian. For a given ao, Eq. (2) is solved 
in the BO-approximation, yielding the field-dressed BO 
curves Pei,i(P, ao) and the dressed energies bFi,j/(ao) [41], 
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with the indices i = 1, 2,... and ^ = 0,1,... denoting the 
electronic and vibrational states. To treat pulsed laser 
fields, we let the maximal quiver amplitude vary with 
the field envelope, Oq Q!o(f) = cko5(f)- The lowest 
BO curve is plotted in Fig. 2 for ao(t) — 0,1.87 and 3.25. 
With increasing ao{t), the BO curve is shifted upwards in 
energy, towards greater R, and becomes gradually shal¬ 
lower. The latter implies that the dressed vibrational 
time scale, T^{ao{t)) = 2 -k j [Wi^i{ao{t)) - lTi,o(ao(f))], 
increases with ao{t) (see also [35]). 

We now present a qualitative model of the dissocia¬ 
tion mechanism. The validity of the model is deter¬ 
mined later by TDSE results. Let r be the time scale 
for the turn-on (and -off) of the laser pulse, and ag.th 
the quiver amplitude satisfying VFi^o(c«o,th) = E^, i.e., 
when the dressed ground state equals the dissociation 
limit [inset of Fig. 1(d)]. Provided the pulse satisfies 

(i) oo > ao,th, (ii) Tv(0) ~ r, and (hi) r < T^iao), the 
dissociation process occurs as follows [Fig. 2(b)]. Dur¬ 
ing the turn-on of the pulse, (ii) ensures the popula¬ 
tion to follow the field-dressed ground state adiabatically. 
At the field maximum, (i) implies that the bound WP 
populates dressed eigenstates with energies greater than 
Ad- Due to (iii), the turn-off of the pulse can be con¬ 
sidered sudden, and the nuclear WP does not feel the 
fast change of the electronic potential, leaving its posi¬ 
tion and energy unchanged. After the pulse, the nuclear 
WP is trapped above E^, resulting in dissociation via 
the field-free electronic ground state, with NKER given 
by EN(ao) = W^i,o(ao) - E^. 

The laser field regime for which Pdis is nonzero in 
Fig. 1(d) satisfies (i)-(iii). For HJ, we have ooyh = 2.41, 
rv(0) = 15.2 fs and T^{ao,th) = 41.7 fs. In Fig. 1(d), 
the onset of dissociation is indeed around the dashed 
line corresponding to ao = ao,th- This agreement sup¬ 
ports the physical picture of the model. The frequency 
range = 2 — 8 in Fig. 1(d) where Pdis is nonzero cor¬ 
responds to r = 9.0 — 17.9 fs, which fulfills condition 

(ii) at least approximately. For og > cto.th, we have 
T <C Tv(ao,th) < Tv(Q!o), and condition (iii) is satis¬ 
fied as well. DI occurs throughout the whole duration 
of the pulse due to higher-order Floquet components in 
V'eN(a; -I- f}a{t),R). In the stabilization regime, DI is 
greatly suppressed, leaving the population trapped above 
Ed to dissociate. This explains Pion + Pdis = 1 in the pa¬ 
rameter regime of Fig. 1(d) where Pdis is nonzero. 

Figures 2(a) and (c) illustrate the cases for which (i) 
is satisfied, but (ii) and (iii) are not. In the case where 
the pulse duration is short, r ^ Tv(0)) the perturba¬ 
tion can be considered sudden, and the initial popula¬ 
tion is unaffected, resulting in no dissociation [Fig. 2(a)]. 
For long pulses r 3> Tv(q!o), the adiabatic approxima¬ 
tion is accurate, and the initial state follows the dressed 
ground state throughout the whole pulse resulting in a 
final bound distribution similar to the initial population 
[Fig. 2(c)]. For the case where (ii) and (iii) are satisfied. 



FIG. 3. (color online).Vibrational distributions at the end of 
the pulse (bottom panels), and NKER spectra for dissociation 
in the ground electronic Isug state of (top panels). The 
laser pulses have lo = 2.278, varying Fq, and different N^- The 
vertical axes are equidistant in energy using the same scale 
for the bound (bottom panels) and continuum (top panels) 
parts. Each row has the same color scale. The dotted lines 
indicate the energies of the field-dressed vibrational ground 
state at field maxima with respect to the dissociation limit, 
Wifi{ao) — Ed- 


but with ao < agyh) no dissociation occurs. Instead, 
an excited vibrational WP is created containing field- 
free vibrational states with quantum numbers v satisfy¬ 
ing Wi^^{0) = IFi,o(ao) [Fig. 2(e)]. 

The mechanisms described in Fig. 2 imply that by 
varying r and og, the vibrational populations, dissoci¬ 
ation yields, and NKER can be controlled. This predic¬ 
tion is confirmed by the TDSE results of Fig. 3, which 
shows the bound vibrational populations and the NKER 
spectra for the Isdg state at the end of the pulse, for 
w = 2.278, Eg = 1 - 17.9, and W = 10 - 500. We 
checked that no dissociation occurs via the excited elec¬ 
tronic states. For W = 10, the i/ = 0 level is most popu¬ 
lated, except at Eg > 14 where strong non-adiabatic cou¬ 
plings excite some vibrations. No dissociation occurs, in 
agreement with Fig. 2(a). For = 140, the vibrational 
population and the NKER spectra follow approximately 
the dotted line corresponding to IFi,g(ag) —Eg. This is in 
agreement with Figs. 2(b), 2(e), and the accompanying 
discussions, where we argued that the energy of the final 
WP with respect to Eg after the pulse is IFi,g(ag) — Eg. 
If IFi,g(ag) < Eg, the WP is bound, and field-free vi¬ 
brational states satisfying ~ IFi,g(Q;o) are popu¬ 

lated. If IFi,g(ac) > Eg, the WP dissociates with NKER 
En ~ Wifi{ao) — Eg. Note that the onset of dissociation 
occurs at Eg ~ w^oggh- From = 140 to iVc = 500, the 
adiabatic approximation for the evolution of the WP de¬ 
scribed in Fig. 2(c) becomes gradually more appropriate, 
with the field-free vibrational ground state populated up 
to increasingly larger Eg. For = 350 — 500, no dissoci¬ 
ation occurs. By suitable choices of r and Eg, we can thus 
control the final vibrational populations and the NKER 
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FIG. 4. (color online). Time-resolved analysis showing the 
buildup of the JES for DI during a laser pulse with parame¬ 
ters ao = 2.41, uj = 2.278, and = 100. Diagonal lines in¬ 
dicate the energy conservations -|- En = Eq+oj (left lines) 
and Ee -I- En = Wifi{ao) + ui. The upper and side subpan¬ 
els show respectively the photoelectron and NKER spectra. 
Qi,..., Q 4 denote the four maxima. 

spectra. 

To understand the interplay between DI and dissocia¬ 
tion in this dynamical regime more deeply. Fig. 4 presents 
the formation of the JES describing the probability of 
measuring electronic Eg and nuclear E-^ energies in the 
DI process, determined as in [42]. Note that compared to 
Fig. 3, the NKER energies from DI are much greater than 
that from dissociation. The JES after the pulse {t = 00 ) 
shows four distinct peaks, denoted by Q^,..., Q 4 . The 
physical picture in Fig. 2(b) explains these. At t = —6.7 
fs, the Stark shift A(t) = Wi^o{ao(t)) — Eq is negligi¬ 
ble, and ionization produces in the JES along the line 
corresponding to the one-photon resonance Eg + Ejq = 
Eq+oj. At t = 0, according to the model, the nuclear WP 
has the largest Stark shift A(0) = 0.09732, and ioniza¬ 
tion leads to Q 2 along the shifted one-photon resonance 
Eg + An = Wi^o{ao) +uj. At t > 0, the population stays 
above Eq and dissociates. During dissociation, ioniza¬ 
tion occurs at large R, leading to Q 3 and Q 4 with low 
En in the JES at t = 6.7. The structure of Q 3 and Q 4 
will not be discussed here. Note that Q 4 is unchanging 
in magnitude and position for t > 0 , consistent with the 
model prediction that no field-free bound states become 
repopulated during pulse turn-off. For this reason the 
dynamic interference effect [43, 44], where ionized WPs 
with equal continuum energy created during the rising 
and falling edges of the pulse interfere, is not observed. 
These insights indicate that dissociation and probing of 
dissociation by DI can be achieved using a single pulse, 
with the nuclear WP being promoted to the field-free dis¬ 
sociation continuum during turn-on (pump) and probed 
through DI during turn-off. 

The dynamics in terms of the physical model hinges 
on two requirements: a) stabilization with respect to DI; 
and b) increase of the Eg\ i and with ag- Consider 
first a). We expect stabilization to occur in two-electron 
molecules such as H 2 , D 2 , and HeH+, provided uj > Ip 


[45, 46]. For other diatomics like N 2 , O 2 , CO, or NO, the 
orbital energies can be grouped into two classes formed 
by innershells with energies less than -11, and outer- 
shells with energies larger than -2. The same classifica¬ 
tion holds for triatomic molecules (e.g., CO 2 and H 2 O), 
and even larger systems like naphthalene. For typical w’s 
considered in this work, e.g. 2-3, the outershells satisfy 
UJ > Ip where we expect stabilization, while innershell 
electrons can be excited to Rydberg states by multipho¬ 
ton absorption and subsequently be stabilized [28, 29]. 
Now consider b). For large og, it has been shown that 
Eg\,i cx Qfg [45, 47], while a perturbative analysis for 
small ag [35] shows that for a general diatomic molecule, 
Eeip, i?g, and Tv increases with ag. For intermediate 
ag, using quantum chemistry codes, the requirement b) 
was also fulfilled for H 2 [47-49]. Even when the molecule 
is not aligned with the linearly polarized field, the men¬ 
tioned two requirements a) and b) still holds [41, 47], in¬ 
dicating that the observed dynamics would also occur in 
non-aligned molecules. In light of these considerations, 
the described dynamics are expected to be observable 
in two-electron diatomics, while further studies are re¬ 
quired to determine whether this is generally the case. 
For a molecule X 2 , with X some atomic species, the 
charged fragments X+ and e from DI can be detected 
in coincidence by cold-target recoil ion momentum spec¬ 
troscopy [50]. The resulting JES will contain structures 
satisfying Eg + 2Ey^+ = kFi,g(ag) -I- oj for outershell DI, 
and Eg + 2Ex+ = kFi,g(ag) -I- riminW for innershell DI, 
with rimin the minimum number of photons required. A 
large spread in Ey^+ containing low-energy X+ fragments 
is a signature of continuous ionization during dissocia¬ 
tion towards larger R, and thus a direct indication of 
the dressed dissociation channel (see Q 2 -Q 4 in Fig. 4 for 
HJ). We believe that an experimental verification would 
be challenging, but possible in the future, involving the 
generation of relevant pulses, molecular alignment, and 
channel-specific measurement techniques. 
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